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ABSTRACT 

We present observations carried out using the 10.4 m Gran Tclescopio Ca¬ 
narias and an interpretative model of the dust environment of activated asteroid 
313P/Gibbs. We discuss three different models relating to different values of the 
dust parameters, i.e, dust loss rate, maximum and minimum sizes of particles, 
power index of the size distribution, and emission pattern. The best model corre¬ 
sponds with an isotropic emission of particles which started on August 1st. The 
size of grains were in the range of 0.1 — 2000 /irn, with velocities for 100 /irn parti¬ 
cles between 0.4 —1.9 m s' 1 , with a dust production rate in the range of 0.2 — 0.8 
kg s' 1 . The dust tails’ brightness and morphology are best interpreted in terms 
of a model of sustained and low dust emission driven by water-ice sublimation, 
spanning since 2014 August 1st, and triggered by a short impulsive event. This 
event produced an emission of small particles of about 0.1 /iin with velocities of 
~4 m s' 1 . From our model we deduce that the activity of this Main-Belt Comet 
continued for, at least, four months, since activation. 

Subject headings: comets: individual (313P/Gibbs) — methods: numerical - 
minor planets, asteroids: general 


1. INTRODUCTION 


The Main-Belt Comet (MBC) 313P/Gibbs, hereafter 313P, was discovered in the course 
of the Catalina Sky Survey’s 0.68-rn Schmidt telescope on September 24.3, 2014. This 
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MBC was desc ribed as a faint object with a narrow tail 20" in length and position angle 
270° ([Gibbs & Sato 2014). These objects have the dynamical properties of asteroids (i.e., 
Tisserand parameter with respect to Jupiter larger than 3), but display comet-like structure 
due to dust and/or gas emission. 

The origin of the activity in these objects is still unclear. A variety of mechanisms 
have been proposed. While some of th em are associated with impulsive events, such as colli- 
sions and rotational break -up (see.e.g-. lMoreno et al.ll20l4 IJewitt et aU2014ai lYang fc Hsieh 


201ll:llshiguro et al. 


201 lh others show an activity sustained over time probably due to water 


ice sublimation (Jewitt et al.l l2014bt iMoreno et al.l 120131: iLicandro et al.l l2013t iHsieh et al. 


2012d ; Moreno et ah 2011a). Of these driving mechanisms of the MBCs activity, it seems 
that the most likely is water-ice sublimation, and an attempt to link this fact with the Earth’s 
oceans has been established. Due to the high temperatures on Earth when it was formed, 
an external source of water is widely accepted. Comets were considered the best candidates 
for this theory, but the measurements of the isotope ratios of deuterium to hydrogen (D/H 
ratio) have shown a wide range of values, but most of them being higher than the oceanic 
D/H ratio (jAltwegg et al.l 120151) . The existence of water in the Main-Belt should not be a 
surprise. Actually, half of the outer belt asteroids show absortion features, which can be 
attributed to the presence of hydrated minerals, and water ice has been detected in (24) 
Thermis. 

N-Body integr ations have shown that MBCs are dynam ic ally stable on time sc ales of 
100 Myr or longer ( Ipatov fc Hahn 1999 ; Hsieh et al. 2012a b; Stevenson et al. 2012), and, 
consequently, they are consi dered native members of the Main Asteroid Belt and not captured 
objects from elsewhere ( Hsieh et al.]2009 ). For fu r ther d etails on those objects we refer the 
readers to the works of Bertini~(l2011 1 and Jewitt ( 2012 ). 


Observations and analysis of 313P have recently been carried out by Jewitt et ah) (120151 ). 
who report an equilibrium sublimation of dirty water-ice from a small patch of the surface as 
origin of the dust tail. The authors ruled out other possible explanations for the origin of the 
activity. Also. lHsieh et al.l (I2015h have studied this ob ject, reaching t he sam e conclusion that 
the activity is likely sublimation-driven. In addition, Hui fe Jewitt (12015 1 presented predis¬ 
covery observations of 313P in 2003 and 2004, establishing that, after 133P/Elst-Pizarro and 
238P/Read, this MBC is the third activated asteroid to show mass loss in different orbits. 
In this work, we report images taken at 10.4 m Gran Telescopio Canarias (GTC) of 313P on 
three different dates. We use our Monte Carlo dust tail code, in order to provide estimates 
of the ejected mass, the particl e ejection velocities , and their size d i stribu tion. Finally, we 
compare our results to those of Hewitt et al.1 (120151) and lHsieh et al.1 (120151)) . 
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OBSERVATIONS AND DATA REDUCTION 


The observations of 313P were carried out on the nights of September 29, November 4 
and December 16 2014 (hereafter 29Sept, 4Nov and 16Dec), using Sloan / and (j filters in the 
Optical System for Image and Low Resolution Integrated Spectroscopy (OSIRIS) camera- 
spectrograph ( Cepa et ah 2000 ; Cepa 2010 1 at the GTC. OSIRIS provides a field of view 
7 '.8 x 7 '.8 and a pixel scale of 0".125 pixel -1 . In order to increase the signal to noise ratio 
the data were binned in 2 x 2 pixels, so that the spatial resolution of the images become 
263 km pixeR 1 , 288 km pixeR 1 , and 373 km pixeR 1 on the observation dates. The images 
were bias and flat corrected using standard techniques and the calibration was performed 
using photometric zero points determined by standard star observations. The r ! images of 
the 313P were transformed into R standard Johnson-Cousins magnitude, where we assume 
the same spectral dependence for both the object and the Sun within the bandpasses of these 
two red Liters. Then, images were converted to solar disk intensity units (SDU) appropriate 
for analysis in terms of the dust tail models. The resulting images and the observational 
details are shown in Fig. |T]and Table [fl respectively. 


THE MODEL 


The analysis was performed using our Monte Carlo dust tail code, which allows us to 
generate synthetic dust tail brightness images to be compared with the observations. This 
code has been su ccessfully used in previous works on characterization of dust e nvironments of 


MBCs (see, e.g., Moreno et al.ll2014].l2013l.l2011a ) and Jupiter Family Comets ( Moreno et al. 


2012b : Pozuelos et al. 2014a 3). In that code, we computed the motion of the particles after 
they had been ejected from the nucleus and submitted to the gravity force of the Sun and 
the radiation pressure, describing a Keplerian orbit around the Sun. The type of orbit of the 
dust particle (elliptica l, hyp erbolic, or parabolic) is a function of its ejection velocity and 


the /3 parameter (e.g., Fulle 1 9891 ). This parameter can be expressed as j3 = C pr Q pr /(2pd ), 
where C pr = 1.19 x 10 -3 kg m -2 , Q pr is the radiation pressure coefficient, p is the particle 


density, and d is the particle diameter. Owing to the many model inputs, we set some of 
them to a specific value. Thus, particles are considered spherical, with density assumed at 
p = 1000 kg m -3 , and their refractive index is set at 1.88 + O.Tli, which corresponds with 


carbonaceous composition (lEdohl 119831 ). Using the Mie theory for spherical particles, we 
find th at geo metric albedo is p v = 0.04, and that Q pr ~ 1 for particles of radius r > 1 pm 
(IMoreno et al.l (j2012at ), their Figure 5). 


The particle size distribution is assumed to follow a time-dependent power law with 
index a. This index is within the range of -4.2 to -3, found in a number of comets (see, e.g., 
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Jockers (11997 ) and the references therein) The particle ejection velocities are described as a 


function of /3 parameter as v(/3) = d 0 /3 7 . This relationship is generally accepted for terminal 
velocities of comet dust, and also for fragments ejected after collision experiments (see, e.g., 
Giblinl 1998; Onose & Fuiiwara 2004!) • Thus, 7 , and the time-dependent parameters, i.e., v 0 , 
the dust mass loss rate, the size of particles (r max and r min ), and the index a of the particle 
size distribution as a function of the heliocentric distance are fitting parameters. 


4. RESULTS 


In order to fit the available observations of 313P, all the parameters described above 
must be introduced as inputs of the code. In order to explore the largest number of possible 
scenarios, we set 7 to three different values, i.e., 7 = 1 / 2 , 1 / 8 , and 1 / 20 , and looked for 
the best fit model in each case, varying the rest of the parameters. The general procedure 
consists of a trial-and-error procedure, starting from the most simple scenario, where we 
consider an isotropic ejection outgassing model, having minimum and maximum particle 
radius of r m) ; n = 1 /im, r max = 1 cm, and set a = —3.5. Regarding the dust mass loss rate 
and Do, we impose a monotonically symmetric evolution with respect to perihelion. Then, 
we start to vary the parameters until an acceptable agreement with the dust tail images is 
reached. 


Our first model corresponds with 7 = 1/2 (model I). This is con sistent with pre¬ 
vious studies of the MBCs P/2010 R2 (La Sagra) ( Moreno et al. 2011a ) and P/2012 T1 
(PANSTARRS) ( Moreno et al. 2Q13). After a number of attempts, the best fit model con¬ 
sists of an emission of particles with sizes ranging between 1-2000 /mi. The power index 
of the size distribution is set to a constant value of a = —3.2. The terminal velocities for 
particles of 100 /mi are between 0.6-5 m s -1 . The maximum dust production rate is 1 kg s -1 
and the total dust emitted is estimated to be 3.2xlO 6 kg. The emission of particles started 
at r h = —2.396 AU, i.e., August 1 st, and the maximum values for the dust production rate, 
terminal velocities and the r max are reached at r h = +2.395 AU, i.e., September 21st. 

The second pro posed model corresponds to 7 = 1/8 (model II). This assumption was 
also used in Moreno et a h (2014) where the authors studied the activated asteroid P/2013 
P5 (PANSTARRS). The range of size particles is 0.1-2000 /mi, with the values of the power 
index of the size distribution between a = —3.1-(—3.3). The terminal velocities are in the 
range of 0.4-1.9 nr s _1 for particles of 100 /mi. The peak of the dust production rate is 0.8 
kg s -1 and the total dust ejected is 3.4xlO 6 kg. In this model, the emission started at the 
same time as the model of 7 = 1 / 2 , i.e, at 77 = —2.396 AU (August 1st) and the maximum 
value of the dust parameters were reached also at r h = +2.395 AU (September 21st). Despite 





















5 


starting and reaching the maximum values for the dust parameters at the same time as the 
model I (August 1st and September 21st respectively), the evolution of those parameters as 
a function of the heliocentric distance were different. The terminal velocities for particles 
of 100 /iin at the beginning of the emission reached the same values as those reached at 
77 = +2.395 AU (September 21st). In addition, the smallest grains (0.1 /jrn) are ejected just 
at the start of the activity. This points to an emission of particles at high velocities at the 
beginning of the activity. 


The last model in our study set 7 = 1/20 (model III). This option has been used 
before in the study of the MBC P/2012 F5 (Gibbs) ( Moreno et al. 2012a ) and (596) Sheila 
( iMoreno et al.l 1201 lbt ). The size of particles were in the range of 0.1-2000 /iin and a took 
values between — 3.1-(—3.2). In this case, the terminal velocities of 100 //rn particles were 
slower than in the other two models, with values of 0.4-0.9 m s' 1 . The maximum dust 
production rate reached the value of 0.80 kg s' 1 , being the total emitted dust 2.9xlO 6 kg. 
As before, the start of the activity seems to be at 77 = —2.396 AU (August 1st), but in this 
case, the maximum values of the dust parameter were reached at perihelion, i.e., August 
28th (77 = ±2.3916 AU). 


I 11 order to determine which of the proposed models give the best fit, we computed the 
fitting quality to the observations by the quantity a = y/^JJobs — Imod) 2 / N, where N is 
the number of pixels of the images, I 0 b s and I mo d are the observed and fitted images. The a 
parameters at each model and date are given in Table [2] In addition, in order to facilitate 
an easy inspection of the models, we display in Figure[2]the comparison between models and 
observation on November 4.05. We see that model II gives the lowest value of a at all three 
dates of observation, and, therefore, this model is favored against models I or III, although 
none of those models can be excluded. I 11 Figure [3] we display the evolution of the dust 
parameters as a function of heliocentric distance, and in Figure [4] we present the comparison 
of this model with the observational data. 


5. DISCUSSION 

From our analysis of the dust environment of 313P, we have shown that the model 
having 7 = 1/8 is favored respect to 7 = 1/2 or 7 = 1/20 models. This indicates a very 
weak dependence of v{/3) with (3. The size of particles is in the range of 0.1-2000 /irn. 
The smallest grains in that range were ejected just at the beginning of the emission at 
rh = —2.396 AU (August 1st) with terminal velocities of 4.4 m s' 1 , i.e., ejection of small 
particles at moderately high velocities at the beginning of the activity. In general terms, the 
peak of the activity is reached at 77 = ±2.395 AU (September 21st). The dust production 
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rate was in the range of 0.2-0.8 kg s -1 , which is weak activity, but it lasted about four months 
from the start of the emission until our last observation on December 16.85, with a total of 
dust emitted of 3.4 x 10 6 kg. This emission is consistent with the sublimation of volatile 
ices found in previous works, such as P/2005 U1 (Read) in 2005, when the comet displayed 
a substantial coma with large particles of r > 100 /mi, and ejection velocities of ~ 0.2 — 3 


-.-l 


m s " OHs i eh et al.l 120091 ). Therefore, we conclude that 313P activity is sustained probably 


by water-ice sublimation triggered by a short impulsive event, like a collisional excavation 
of surface ice. 


Jewitt et ah (2015), obtain a continuous mass-loss at rates of ~0.2-0.4 kg s -1 , with size 
particles of 50-100 /irn. The maximum values of those ranges are a bit lower than ours, but 
close to them. The authors concluded that the ejection of material is consistent with the 
action of sublimation of dirty water-ice from an exposed patch covering only a few hundred 
square meters, and protracted from September until November. In our model we do not 
constrain the activity to an active area, but the emi s sion of particles can be explained by 
an isotropic ejection model. In addition, Hsieh et ah (12015 ) gave a total of dust ejected in 
the range of 2-4x10' kg, which is larger than the dust production obtained in our model. 


Regarding the origin of the activity of the MBC, both Jewitt et ajJ (120151 1 and Hsieh et ah 


(120151 ). suggest an impact or local avalanche as the nature of the trigger to expose ice in first 
place, which is consistent with the results obtained in our model. 


In addition, 1 Ini & Jewitt (120151 ) report the existence of archival data showing that 
313P was active in its perihelion passage in 2003, favoring the theory of sublimation of 
water-ice. This puts 313P in the same category of active asteroids as 133P/Elst-Pi zarro 


and 2 38P/Read, which exhibit activ ity at more than one perihelion passage (IHsieh et ah 


20041 . 1201 ll . 120131: iMoreno et al.l 120131 ) . This fact is considered the most reliable indicator of 


water-ice sublimation. Despite this evidence, until today, no gas has been spectroscopically 

ration incomplete (see, e.g., 
Therefore, the driver of the 


de Val-Borro et al. 

_--2—, 

2012 ; 

Licandro et al. 

2011 ; 

Jewitt 

2009) 


activity in 313P can not be clearly established. 


6. CONCLUSIONS 

The application of the Monte Carlo dust tail model to images of the recently discovered 
MBC 313P/Gibbs, carried out at GTC, has revealed that the observed morphology and 
evolution of the tail is consistent with the emission of particles started on August 1st (r^ = 
—2.396 AU), with particle sizes in the range of 0.1-2000 /mi, which were ejected isotropically. 
The dust loss rate varies between 0.2-0.8 kg s -1 , reaching the peak of the activity at r h = 
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+2.395 AU (September 21st). The terminal velocities are parametrized by v{j3) = no/3 7 , 
where we favor 7 = 1/8 over the other tested values (7 = 1/2 and 7 = 1/20). Those 
velocities showed that at the beginning of the emission, the smallest grains (0.1 /mi) were 
moderately rapidly dispersed, with velocities of 4.4 m s -1 . The total dust ejected is found 
to be 3.4xlO 6 kg. Therefore, we conclude that the sustained low activity observed for a few 
months is likely driven by water-ice sublimation, due to a collisional event, which exposed 
fresh ices to the solar radiation. This scenario and our estimates of the dust production rates, 
ejection velo c ities, and particle sizes are in agreement with the results recently reported by 


Jewitt et al.l ( 20151 ) and iHsieh et al.l ( 2015 ) for this MBC. 
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Fig. 1.— Panel (a) displays the 313P/Gibbs imaged on 2014 September 29.07. Isophote 
levels (solid blue lines) are 0.25 xlO -14 , 0.65 xl0~ 14 , 1.50 xlO” 14 , 0.35 xlO -13 , and 1.50 xlO -13 
(SDU). The scale is 263 km pixel -1 . Panel (b) shows the same as (a) but for the 2014 
November 4.05. Isophote levels are 1.00 xlO -14 , 1.50 xlO -14 , 3.00 xlO -14 , and 6.00 xlO -14 
(SDU). The scale is 288 km pixel -1 . Panel (c) corresponds with 2014 December 16.85. 
Isophote levels are 0.28xl0 -14 , 0.5xl0 -14 , l.OOxlO -14 , and 1.75xl0 -14 (SDU). The scale is 
373 km pixel -1 . In all images North is up and East is to the left. 
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Fig. 2.— Comparison between models on the observation date of 2014 November 4.05. 
Panels (a.l) and (a.2) correspond to the model with 7 = 1/2, (b.l) and (b.2) with 7 = 1/8, 
and (c.l) and (c.2) with 7 = 1/20. In all cases, the black line is the observed image and the 
red line the fit. The blue line along the tails in panels (a.l), (b.l) and (c.l) correspond to 
the brightness scan displayed in (a.2), (b.2) and (c.2). The isophote levels are 1.00xl0~ 14 , 
1.50xlO -14 , 3.00xl0“ 14 , and 6.00xlO -14 (SDU). The vertical bars correspond to 3000 km 
in the sky. 
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Fig. 3.— Evolution of the dust parameters obtained in the model II (7 = 1/8) versus the 
heliocentric distance for MBC 313P/Gibbs. In panel a) dust mass loss rate [kg s' 1 ]; b) 
terminal velocities for 100 /mi glassy carbon spheres [m s' 1 ]; c) maximum and minimum size 
of ejected particles [/mi]; d) power index of the size distribution a. The dashed vertical line 
corresponds to perihelion. 
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Fig. 4. — Isophote field and brightness scans comparison between observations and model 
II. The black lines correspond to the observations and the red ones to the model. The 
bine lines along the tails correspond to the brightness scans. Panel (a.l) and (a.2) re¬ 
fer to image on September 29.07 2014. The isophote levels are 0.25 xlO -14 , 0.65 xl0~ 14 , 
1.50xl0 -14 , 0.35xl0~ 13 , and 1.50xl0~ 13 (SDU). Panel (b.l) and (b.2) correspond to the ob¬ 
servation on November 4.05 2014. Isophote levels are 1.00 xlO -14 , 1.50xl0~ 14 , 3.00xlO” 14 , 
and 6.00xlO -14 (SDU). Panel (c.l) and (c.2) refer to the observation on 2014 December 
16.85. Isophote levels are 0.28xl0~ 14 , 0.5xl0~ 14 , l.OOxlO" 14 , and 1.75xl0" 14 (SDU). 
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Table 1. Log of the Observations. 


Observation Date 
(UT) 

Days from 
perihelion 3. 

Th 

(AU) 

A 

(AU) 

Phase 
angle (°) 

Position 
angle (°) 

Resolution 
(km pixeR 1 ) 

2014 Sept. 29.07 

31.6 

2.398 

1.429 

7.9 

317.7 

263 

2014 Nov. 4.05 

67.5 

2.419 

1.565 

14.8 

43.5 

288 

2014 Dec. 16.85 

110.3 

2.464 

2.030 

22.7 

63.2 

373 


a All these observations are post-perihelion. 


Table 2. Fitting quality of the observed images. 


Model 

7 

2014 Sept. 29.07 

2014 Nov. 4.05 

2014 Dec. 16.85 

I 

1/2 

1 .2xl0~ 14 

9.7xl0 -15 

3.3xl0” 15 

II 

1/8 

8 .2xl0” 15 

8 .1xl0” 15 

2.3xl0~ 15 

III 

1/20 

2 .2xl0 -14 

1 .2xl0 -14 

2 .6xl0~ 15 







